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ABSTRACT 

We consider a sample of ten GRBs with long lasting (> 10 2 sec) emission detected 
by Fermi/LAT and for which X-ray data around 1 day are also available. We assume 
that both the X-rays and the GeV emission are produced by electrons accelerated at 
the external forward shock, and show that the X-ray and the GeV fluxes lead to very 
different estimates of the initial kinetic energy of the blast wave. The energy estimated 
from GeV is on average ~ 50 times larger than the one estimated from X-rays. We 
model the data (accounting also for optical detections around 1 day, if available) to 
unveil the reason for this discrepancy and find that good modelling within the forward 
shock model is always possible and leads to two possibilities: either the X-ray emitting 
electrons (unlike the GeV emitting electrons) are in the slow cooling regime or ii) the 
X-ray synchrotron flux is strongly suppressed by Compton cooling, whereas, due to the 
Klein-Nishina suppression, this effect is much smaller at GeV energies. In both cases 
the X-ray flux is no longer a robust proxy for the blast wave kinetic energy. On average, 
both cases require weak magnetic fields (10 -6 < eb < 10 -3 ) and relatively large 
isotropic kinetic blast wave energies 10 53 erg < E 0; kin < 10 55 erg corresponding to large 
lower limits on the collimated energies, in the range 10 52 erg < E^kin < 5 x 10 52 erg 
for an ISM environment with n ~ 1cm -3 and 10 52 erg < E^kin < 10 53 erg for a wind 
environment with A* ~ 1. These energies are larger than those estimated from the X- 
ray flux alone, and imply smaller inferred values of the prompt efficiency mechanism, 
reducing the efficiency requirements on the still uncertain mechanism responsible for 
prompt emission. 
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1 INTRODUCTION 

Gamma-Ray Bursts (GRBs) have two phases. A short 
prompt phase in which the emission is predominantly in 
the sub-MeV range is followed by a multi-wavelength af¬ 
terglow that can be observed up to several years after 
the prompt emission. According to the generally accepted 
internal-external model (Sari V Piran 1997) the prompt 
emission is thought to originate from internal energy dis¬ 
sipation that takes place within the relativistic outflow. 
Conversely, the afterglow radiation is thought to originate 
from interactions between the outflow and the external 
medium (Rees V Meszaros 1992; Paczynski and Rhoads 
1993; Meszaros et al. 1997; Sari V Piran 1997). While the 
origin of the afterglow radiation has been understood in its 
general properties of synchrotron emission (Sari et al. 1998), 
the nature of the prompt radiation has not yet been com- 
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pletely unveiled. Indeed, the nature of both the dissipation 
and the radiative mechanisms is still uncertain. The com¬ 
position of the outflow, its initial energy content, and the 
processes at the origin of the prompt radiation are among 
the most relevant open issues in GRB studies. 

A critical issue in this model is the ratio of energy in¬ 
volved in the first (prompt) and second (afterglow) phases. 
This ratio reflects the efficiency of the prompt radiation pro¬ 
cess, a very important clue on the nature of this mechanism. 
The kinetic energy left in the blast wave after the prompt 
phase can be estimated from afterglow observations. Assum¬ 
ing that the cooling frequency lies below the X-ray band, 
the X-ray luminosity provides a robust estimate of the en¬ 
ergy stored in the accelerated electrons, which in turn is 
directly related to the kinetic energy of the outflow (Kumar 
2000; Freedman V Waxman 2001). Under this assumption, 
several studies have exploited both pre-Swift X-ray obser¬ 
vations (Berger et al. 2003; Lloyd-Ronning V Zhang 2004), 
and Swift X-ray observations (Berger 2007; Nysewander et 
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al. 2009). Most of these studies have inferred a relatively low 
kinetic energy, which implies quite a large prompt efficiency: 
e 7 > 0.5. The discovery of the X-ray plateaus in many of the 
Swift GRBs increased the severity of the efficiency problem. 
The X-ray flux at the beginning of the plateau phase (around 
500 sec) is lower by a factor ~ 3 as compared with the flux 
estimated by extrapolating backwards in time the observa¬ 
tions at ~ 1 day and therefore leads to an estimate of the 
kinetic energy lower by the same factor and to efficiencies of 
up to 90% (Granot et al. 2006; Ioka et al. 2006; Nousek et 
al. 2006; Zhang et al. 2007). 

Internal shocks are the dominant internal dissipation 
process for matter dominated flows (Narayan et al. 1992; 
Rees V Meszaros 1994). Since their efficiency is rather low 
(Kobayashi et al. 1997; Daigne V Mochkovitch 1998; Be¬ 
loborodov 2000; Kobayashi V Sari 2001; Guetta et al. 2001) 
it was expected that after the prompt phase most of the en¬ 
ergy would remain as bulk kinetic energy of the blast wave. 
Alternative scenarios, such as magnetic reconnection (Usov 
1992) (that require a magnetically dominated rather than 
a matter dominated outflow) may reach higher efficiencies, 
leaving less energy in the blast wave. Thus the high effi¬ 
ciency, implied by the X-ray afterglow observations, is gen¬ 
erally considered as a major problem for the internal shocks 
model and suggested that other mechanisms, such as mag¬ 
netic reconnection take place. 

However, two implicit assumptions have been made 
when deriving these estimates: first it was assumed that the 
electrons emitting at X-rays are fast cooling (i.e. the X-ray 
is above all the typical synchrotron break frequencies) and 
second the X-ray flux of these electrons is not suppressed by 
Synchrotron-Self Compton (SSC) losses (Sari V Esin 2001; 
Fan V Piran 2006). If either one of these assumptions is 
not satisfied the energy estimates based on the X-ray fluxes 
might be wrong. 

Observations at higher energies could be helpful in con¬ 
straining the location of the synchrotron cooling frequency 
and assess the importance of the SSC mechanism. GRBs can 
now be studied at energies between 0.1 GeV and 300 GeV 
thanks to the Large Area Telescope (LAT) on board Fermi. 
While only a small fraction of Fermi detected GRBs have 
also been detected by LAT, some of these cases (about 10 
events) are particularly interesting, since they show charac¬ 
teristics suggestive of an external shock origin for the GeV 
radiation: first, the onset of these emission is delayed rela¬ 
tive to the onset of the prompt sub-MeV emission (Abdo 
et al. 2009a); second the LAT component extends long af¬ 
ter the prompt sub-MeV emission stops and third the flux 
of this long lasting component decays as a power-law in 
time. Indeed, these observations are compatible with expec¬ 
tations from forward shock radiation (Kumar & Barniol Du¬ 
ran 2009, 2010; Ghisellini et al. 2010; Wang et al. 2013; Nava 
et al. 2014). 

We can exploit the afterglow observations at high en¬ 
ergies to address the questions of the determination of the 
synchrotron cooling frequency and the importance of the 
SSC mechanism (which leads to a determination of the en¬ 
ergy content of the fireball and the efficiency of the prompt 
mechanism). We examine a sample of GRBs detected both by 
LAT and XRT, and use both the GeV and the X-ray fluxes 
to estimate the blast wave kinetic energy. We show that the 
two estimates are inconsistent with each other. The energy 


inferred from the GeV emission is much larger than that es¬ 
timated from the X-rays. This can be explained if either i) 
the X-ray emitting electrons, unlike those emitting at GeV 
energies, are in the slow cooling regime, or ii) if electrons 
radiating at X-rays are significantly cooled by SSC, while 
those radiating at GeV energies are in the Klein-Nishina 
(KN) regime and cool only by synchrotron. In both scenar¬ 
ios the X-ray flux is no longer a proxy for the blast wave 
kinetic energy. We examine afterglow models for the GeV, 
X-ray and (when available) optical data in order to deter¬ 
mine if one or both of the mentioned scenario represents a 
viable solution to the XRT/LAT apparent inconsistency. 

The nature of the solution depends on the cooling fre¬ 
quency and on the value of the Compton parameter. Both 
depend strongly on the fraction of energy stored in the down¬ 
stream magnetic field cb and on the density of the external 
medium n. Modelling of the data allows us, therefore, to con¬ 
strain these parameters. Recent works on GRBs with GeV 
emission (Kumar V Barniol Duran 2009, 2010; Lemoine et 
al. 2013) and without it (Barniol Duran 2014; Santana et 
al. 2014; Zhang et al. 2015; Wang et al. 2015) suggested 
a distribution for the es parameter that extends down to 
very small values: 10 -6 — 10 -7 . In this work we derive upper 
limits on cb as a function of the external medium density 
and find similar results: upper limits on derived form our 
sample range from 10“ 2 down to 10 -6 . 

The paper is organized as follows. In §2 we describe the 
sample of GRBs used in this paper. In §3 we show that the 
kinetic energies of the blast waves inferred from GeV data 
are much higher than those inferred from the X-ray obser¬ 
vations, and in turn the corresponding prompt efficiencies 
are much smaller. We examine, in §4, both scenarios and 
we provide analytic estimates of the parameters needed to 
explain the observations. We present the numerical method 
used to model the data and describe the results of the de¬ 
tailed numerical modeling in §5. Finally, in §6 we summarize 
our conclusions. 


2 THE SAMPLE 

We consider a sample of GRBs detected both by Fermi-LAT 
and by Swift-XRT. We included in our sample only those 
bursts for which the LAT emission lasted much longer than 
the prompt phase, since we are interested in cases in which 
the LAT is most likely afterglow radiation from the external 
forward shock. When available we use optical observations 
to further constrain the solution. The final sample includes 
ten GRBs, listed in Table 1. The redshift has been firmly 
measured for nine of those. 

To ensure that the measured fluxes are most likely orig¬ 
inated at the forward shock with no reverse shock contribu¬ 
tion and to make a proper comparison with previous studies, 
we chose XRT observations that are as close as possible to 
~ 1 day. We select two observation times for each burst. The 
first such that it satisfies the following conditions: i) they are 
as close as possible to one day, ii) subsequent to the end of 
any plateau phase in the light-curve that might be present 
(this phase is present only for one burst in the sample, GRB 
090510) and iii) before the jet break time (when present). 
The second observation time is chosen such that it is as far 
removed from the first observation as possible (in order to 
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enable a good estimate of the light-curve power law index) 
while satisfying conditions ii) and iii). Such late observa¬ 
tions are not available for the LAT and we therefore use 
earlier epochs for the LAT measurements. For LAT we use 
two data points as far removed in time as possible given that 
they are after Tqo (the duration of the prompt emission as 
determined by the lower energy Fermi / GBM detector) and 
provided that the flux measurement is still well constrained. 

For the LAT fluxes we used values reported in the First 
Fermi/LAT GRB catalog Ackermann et al. (2013a) for nine 
of the bursts, while for GRB 130427A we used the results 
reported in Ackermann et al. (2014). The X-ray fluxes are 
taken from the Swift/XRT GRB light curve repository 1 
(Evans et al. 2007, 2009) and the optical fluxes are collected 
from the literature. The complete list of LAT, X-ray and 
Optical fluxes and observations times is reported in Table 1. 


3 KINETIC ENERGY OF THE BLAST WAVE 

We begin by re-visiting the estimates of the kinetic energy 
of the afterglow that are based on X-ray observations. As 
in previous works, we assume for now that the X-ray ra¬ 
diation around one day is above the cooling frequency v c 
(the frequency above which electrons cool efficiently) and 
that Compton losses are negligible (i.e. 7 C 1, where Y 
is the Compton parameter). In this case the kinetic energy 

4 2 — p 4(1 — p) 3p — 2 

scales as: E^ k oc F^- +p e^ +p e e 2+p t£ +P , where tx is the 
time of X-ray observations and Fx is the X-ray flux at that 
time (see Eq. 1 below). To estimate the kinetic energies, we 
use the data in Table 1 (in particular we use tx,l and Fx, i 
as well as tGeV ,i and Fo e v, i) and assume that e e — 0.1, 
p = 2.5 and cb = 10 -2 . The results depend only weakly 
on p and cb■ The resulting energy estimates are similar to 
those found in previous studies (Frail et al. 2001; Panaitescu 
and Kumar 2001a,b; Granot et al. 2006; Ioka et al. 2006; 
Nousek et al. 2006; Zhang et al. 2007; Nysewander et al. 
2009; Cenko et al. 2011). In particular, they imply large 
gamma-ray prompt emission efficiencies: (e 7j x) = 0.87 (see 
Table 2), where e 7 = E 1 /(E 1 + E 0j kin), and the sub-index 
indicates whether X-ray (“X”) or GeV fluxes (“GeV”) are 
used to estimate it. The efficiency is as large as 0.99 for 
GRB090902B. 

We repeat the same calculation using the GeV fluxes, 
assuming that this high energy radiation is also produced 
by forward shock synchrotron. We use photons detected by 
LAT in the 0.1-10 GeV range and it is therefore expected 
that they are in the fast cooling regime of the spectrum. If 
this is the case and if e e and the prompt efficiency are nar¬ 
rowly distributed, then a strong correlation between GeV 
luminosity and E 1 is expected (Ghisellini et al. 2010; Beni- 
amini et al. 2011). This has indeed been found by Nava et 
al. (2014) for a sample of 10 GRBs which almost completely 
overlap with the sample of bursts studied in the present 
work, strengthening the hypothesis that the high energy ra¬ 
diation is synchrotron radiation from the external forward 
shock at frequencies larger than v c . Compton cooling of the 
emitting electrons is suppressed by the KN lower cross sec¬ 
tion (Wang et al. 2010), that is relevant at such high en- 

1 http://www.swift.ac.uk/xrt_curves/ 


ergies. Therefore, the GeV flux is expected to be a better 
proxy of the kinetic energy than the X-ray flux. 

The kinetic energies Eq^ u obtained from the GeV flux 
are much larger than kinetic energies EQ kin obtained from 
X-rays (see Table. 2). The inconsistency between the esti¬ 
mated kinetic energies can be better appreciated in terms 
of their ratios, shown in Fig. 1. Typically Eq^ u is at least 
ten times and up to a hundred times larger than E^ kin . 
The efficiencies inferred using the GeV observations are con¬ 
siderably lower than those obtained from the X-ray data: 
(e 7 ,Gev) = 0.14, as compared to (e 7 ,x) = 0.87 inferred from 
X-ray data. 

This ratio, FQ kin /FQjY n ^ i s independent of the as¬ 
sumptions made on cs,e e (assuming these quantities do 
not change significantly between the time of the two ob¬ 
servations), and it depends very weakly on p through: 

3p — 2 2 p 

( tx/tGeV ) 2+p {yx/vGev) 2+p (where vx — 1 keV is the fre¬ 
quency at which the X-ray flux is obtained and VGeV — 
100MeV is the equivalent frequency for GeV). We estimate 
the ratio Eg kin /E$fV n for p — 2.5 and incorporate the un¬ 
certainty on p (2.1 < p < 2.8) in the error bars in the figure. 
The error bars also account for possible radiative losses that 
decrease the kinetic energy content of the fireball between 
the time of the GeV and the X-ray observations. To account 
for these losses, we follow Nava et al. (2013) to estimate 
the radiative losses for each burst assuming e e — 0.1 and 
a ’’fast cooling” regime (this is a conservative assumption: 
it assumes that all the energy dissipated at the shock and 
transferred to the electrons, is radiated). It is clear from Fig. 
1 that the variation of p and possible radiative losses are not 
sufficient to explain the difference. The inconsistency implies 
that at least one of the naive assumptions was wrong. We 
turn in the next section to modeling of the afterglow that 
will enable us to examine this issue. 


4 ANALYTIC ESTIMATES 

We turn now to model the GeV and X-ray data in order 
to understand which one of the initial assumptions is not 
valid. To recall, these two assumptions were that either both 
the X-rays and the GeV emission are in the fast cooling 
regime and that neither emitting electrons are cooled by 
Inverse Compton. We envisage two possibilities. In the first, 
“SSC suppressed ” scenario, both the X-rays and the GeV 
photons are above v c , but the X-ray flux is suppressed due 
to IC losses, while the GeV emitting electrons are in the KN 
regime and their flux is not suppressed (Wang et al. 2010). 
If this is the case, and if Compton cooling is not accounted 
for properly, then a small kinetic energy is inferred from the 
relatively low X-ray flux, while a higher kinetic energy is 
inferred from the high-frequency, unsuppressed, GeV flux. A 
second possibility is the “slow cooling” scenario in which the 
X-ray band is below v c . In this case the X-ray flux depends 
strongly also on es and on the external density and much 
higher energy is needed to produce a given flux than if the 
X-rays were above v c . In both scenarios (SSC suppressed or 
slow cooling) the X-ray flux is no longer a good proxy for 
the kinetic energy of the blast wave. 

We begin by considering simplified analytic estimates to 
explore if and how the two proposed scenarios can explain 
the apparent kinetic energy inconsistency. The analytic es- 
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Burst 

t GeV, 1 

F GeV, 1 

t GeV, 2 

F GeV, 2 

t X,l 

F X,1 

l X,2 

F X ,2 

topt , 1 

F opt , 1 

^opt, 2 

F opt, 2 

Z 

t jet 

ref. 


10 — ^ days 

nJy 

10 — Mays 

nJy 

days 

nJy 

days 

nJy 

days 

pjy 

days 

pjy 


days 


080916C 

4.6 

4.9 

1 

72 

6.94 

9.7 

11.57 

3.6 

1.39 

5.5 

3.47 

1.5 

4.35 

> 15.3 

1 

090323 

4.3 

14 

1.6 

54 

2.495 

23 

5.78 

6.3 

1.85 

14 

5.1 

2.7 

3.57 

> 10 

2 

090328 

14 

1.2 

2 

8.1 

3.47 

20 

6.9 

11 

1.63 

25 

2.6 

11 

0.73 

> 10 

2 

090510 

1.3 

16 

0.1 

210 

0.14 

22 

0.062 

100 

1.16 

1.8 

0.14 

9 

0.9 

> 0.75 

3 

090902B 

7.6 

15 

0.2 

400 

0.928 

21 

2.38 

3.7 

1.43 

10 

2.52 

5.7 

1.822 

20 

4 

090926A 

3 

15 

0.4 

42 

2 

80 

11.57 

5.4 

3 

37 

6.1 

14000 

2.1 

10 

2 

100414 

3.2 

7.1 

0.5 

66 

2.3 

7.6 

7.17 

0.33 

- 

- 

- 

- 

1.37 

> 7.4 

- 

110625A 

5 

38 

3 

42 

0.46 

130 

0.139 

1300 

- 

- 

- 

- 

* 

> 0.47 

- 

110731A 

3.5 

0.75 

0.2 

350 

1.22 

61 

3.84 

12 

0.026 

70 

0.012 

240 

2.83 

> 7.5 

5 

130427A 

63 

2 

3 

110 

1.19 

2600 

11.54 

100 

1.15 

130 

0.22 

1000 

0.34 

> 180 

6 


Table 1 . The GeV, X-ray and optical fluxes and observation times for bursts in the sample. The GeV data is taken from the first 
Fermi/LAT GRB catalog (Ackermann et al. 2013a), apart from 130427A, which is taken from (Ackermann et al. 2014). The X-ray data 
is taken from the Swift/XRT GRB light curve repository. The jet break estimates for 090902B and 090926A are taken from ? The last 
column lists the references for the optical data: 1 Greiner et al. 2009; 2 Cenko et al. 2011; 3 Nicuesa Guelbenzu et al. 2012; 4 Pandey et al. 
2010; 5 Lemoine et al. 2013; 6 Vestrand et al. 2014. *For GRB 110625A there is no measured redshift and we use a typical value of z = 1. 




Figure 1. Left panel: ratio between the kinetic energy of the blast wave as estimated from X-rays (E* ki ) and the same quantity as 
estimated from GeV flux Both X-rays and GeV are assumed to be above u c at the time of measurement. Note that the ratio 

of the estimates is independent of e#,e e and n. The solid line marks equal energy, the dashed line marks E* ki = O.IEq^ and the 
dot-dashed line marks E x kin — 0-01 £& Errors on the ratio take into account possible radiative losses (that might reduce the energy 
content of the blast wave between the time of GeV observations and the time of X-ray observations) and different values of the electron 
distribution power law index p (from 2.1 to 2.8). Right panel: efficiencies of the prompt phase of the GRBs in the sample assuming the 
kinetic energies from X-rays (dashed lines) and GeV (solid lines) shown in the left panel. 


Burst 

^ 7,54 

E x 

^ 0 ,^ 71,54 

F GeV 
^ 0 , kin, 54 

e i,x 

£ 7 ,GeV 

080916C 

3.48 

0.6 (1.38) 

15.2 (20.1) 

0.85 (0.71) 

0.19 (0.14) 

090323 

3.44 

0.28 (0.61) 

27.4 (37-3) 

0.92 (0.85) 

0.11 (0-08) 

090328 

0.1 

0.03 (0.07) 

1.08 (1.71) 

0.77 (0.58) 

0.08 (0-05) 

090510 

0.04 

0.001 (0.002) 

0.88 (1.12) 

0.98 (0.96) 

0.04 (0-03) 

090902B 

2.53 

0.03 (0.06) 

21.4 (32-1) 

0.99 (0.98) 

0.1 (0.07) 

090926A 

1.75 

0.3 (0.67) 

9 (12.3) 

0.85 (0.72) 

0.16 (0.12) 

100414 

0.49 

0.023 (0.05) 

2.3 (3.2) 

0.95 (0.9) 

0.17 (0-13) 

110625A 

0.18 

0.02 (0.05) 

11.4 (16.6) 

0.88 (0.78) 

0.015 (0.01) 

110731A 

0.49 

0.2 (0.4) 

1.1 (1.5) 

0.7 (0.54) 

0.3 (0.24) 

130427A 

0.8 

0.15 (0.35) 

2.5 (4.6) 

0.83 (0.7) 

0.24 (0.15) 


Table 2. Comparison of prompt gamma-ray energy (Ackermann et al. 2013a) and derived kinetic energies using either X-ray or GeV 
data for the bursts in the sample. Column 2 lists the prompt isotropic equivalent gamma-ray energies of the bursts. Columns 3,4 list the 
derived kinetic energies using the X-ray flux - at tx ,l (see Table 1) and GeV flux - at tGeV ,l (see Table 1) accordingly. Columns 5,6 list 
the derived efficiencies for the same data. These values are estimated for p = 2.5 and cb = 10 —2 in an ISM medium. They change only 
mildly for different values of p and cb and for a wind medium. The values in parenthesis correspond to the same quantities in case the 
outflow was affected by radiative losses. 
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timates allow us to understand what are the conditions that 
the free parameters (and in particular cb and n) must sat¬ 
isfy in order to reconcile the kinetic energy estimates from 
both X-ray and the GeV observations. 

In order to present analytic estimates we make the sim¬ 
plified assumption that the GeV radiation is always above 
v c and that due to KN suppression, Compton cooling does 
not affect the spectrum at these frequencies (Wang et al. 
2010), i.e. VGeV > vkn, where vkn is the KN frequency, 
above which seed photons undergo KN suppressed IC scat¬ 
terings, and is given by vkn ~ Tm e c 2 /(h(l + z)j e ) ^ 
rm e c 2 /(/i(l+ ^)min( 7 m , 7 c)) (where h is Planck’s constant, 
r is the bulk Lorentz factor, z is the cosmological redshift 
and 7e is the typical Lorentz factor of electrons in the co¬ 
moving frame which is always equal to or larger than the 
minimum Lorentz factor - either for slow cooling or y c 
for fast cooling). If < y c , the upper limit to the KN 
frequency is constant in time and depends only on e e . For 
e e = 0.1, its value is around ~ 10 18 Hz, well below the LAT 
band. At early times might be larger than y c , and vkn 
is larger that what estimated in the previous case, but still 
below the LAT range. Therefore, in first approximation, the 
GeV flux is a proxy of the kinetic energy. E 0 ,kin can be de¬ 
rived from Eq. 1 (see below), using Y = Yc e v 1 (YceV is 
the Compton parameter for GeV emitting electrons). This 
assumption greatly simplifies the analytical expressions. The 
validity of this assumption is revisited in §5, where we nu¬ 
merically estimate YbeV, instead of assuming YoeV 1. We 
find that in any case its value is of order unity or smaller, 
and the approximation Yc e v 1 that we use in this section 
is good enough to derive order of magnitude estimates. We 
fix the fraction of shocked energy in electrons to e e w 0.1. We 
do this for four reasons. First, contrary to the values of 
e e from afterglow modeling do not vary over orders of magni¬ 
tude and are often consistent with 0.1 (Santana et al. 2014). 
Second, this value is consistent with numerical simulations of 
shock acceleration (Sironi V Spitkovsky 2011). Third, Nava 
et al. (2014) have shown (for a sample of 10 bursts that al¬ 
most completely overlaps with the one used here) that e e 
must be narrowly distributed and that GeV light curves are 
consistent with a value 0.1. Fourth, since the radiated en¬ 
ergy is limited to the energy in the electrons {e e EQ^in), the 
inferred energy from the observed flux is larger for smaller 
values of e e . Therefore it cannot be much lower than 0.1 
in order to avoid un-physically large energy requirements. 
These considerations significantly limit the allowed range of 
this parameter and allow us to fix it without much loss of 
generality. We are left with three free parameters: es, p, and 
n (or A* = A/(5 x 10 n g/cm), where A is the wind param¬ 
eter defined by n = ( A/m p )r~ 2 for a wind environment). 
We may therefore use the expression for the X-ray flux to 
obtain cb as a function of n (or A*), and p. 


4.1 SSC suppressed case 

If Compton losses are important, the flux above the cool¬ 
ing frequency scales as F v (y > v c ) oc Fq 2 ^ )/4 (1 + V) -1 . 

Therefore oc ( 4^7 ) 4/(2+p) ( Y x is the Compton 

parameter for the X-ray emitting electrons). The difference 
in energies can be explained as a difference in the respective 
Compton Y parameters. In the SSC suppressed case, the ex¬ 


pression for the X-ray flux, Fx , is that of synchrotron above 
the cooling frequency (Granot V Sari 2002): 


2-Ur) P — 2 

0.855(p — 0.98)e 1,95p (1 I ;) 4 e e p_1 e B 4 

2-j-p 2 — 3 p 

E ki 4 n , 52 t d £ s dp 28 (^)- p / 2 (l+y^mJy, for ISM 

P~2 

0.0381 (7.11 — p) e 2 ' 76p (1 +z)t e/ -1 e B 4 

2+p 2—3p 

. E kin ,52 *day S d l28 i^~ P 2 for wind 

( 1 ) 


where e e = e 6 |pf, cIl is the luminosity distance, tdays is the 
time since explosion in days, both t and z+l 4 are in the ob¬ 
server frame and we use the notation: q x = q/10 x in c.g.s. 
units here and elsewhere in the text (the numerical coef¬ 
ficients of the ISM and wind cases in Eq. 1 are not very 
different). We assume that X-ray emitting electrons are in 
the Thomson regime, Vm < v c and 2 < p < 3. In this regime, 
the Compton parameter Yx is given by (see Appendix A for 
a derivation): 


Y x 


cb( 3 — p)( 1 + Y x ) 



( 2 ) 


Note that this is an implicit equation, as v c also depends on 
Yx (see Eq. 4 below). The typical frequencies, z/ m , v c (where 
Vyn is the synchrotron frequency of the typical electron en¬ 
ergy) are given by: 


l^m 


and 


3.73 (p — 0.67) 10 15 (1 + 7 /2 El(l 52 

fe 4 /2 Cl' 112 ’ for ISM 

( 3 ) 

4.02 (p - 0.69) 10 15 (1 + 7 /2 El(f 52 

4 2 *da/ S 2Hz ’ fOT Wind 


6.37 (p - 0.46) 10 13 e~ 116p (1 + z)~ 1/2 e ~ 3/2 n d 1 

for ISM 

4.40 (3.45 — p) 10 10 e°' 45p (1 + z)~ 3/2 e B 3/2 A~ 2 

. Elin ,52 ^days + Y x)~ 2 Hz, for wind 


( 4 ) 

Since v c depends on the external density, so does Yx, and 
also Fx- Equating an observed X-ray flux to the flux given 
by Eq. 1 we obtain es(n): 


£b 


4x 10 -5 /(p) (oi) 


2p z -14p-|-12 

p 2 +2p—16 


2(12—p 2 ) 


2 (p—4) (p+2) 


7yip 2 +2p—16 rp p' z +2p— 16 

r GeV -6 r X-B 


(2—p) (p+2) (3p—2) (12—p 2 ) (3p 2 —12p-|-4) (p-|-2) 8(p-2) 

t( 2(£+2p-m 2(p^+2p—16) for IgM 


P-2 

3!~P P _ — 

F 3 ~p 


4 —p 

F 3 ~p 

1 v 


2X10 f(p) ^o e i j ± GeV,-6 ± X,-5 

p—2 p —4 12 — 5p 

77-3 7 <3_p> for Wind 


( 5 ) 

where tcev ,tx are the times of the GeV and the X-ray ob¬ 
servations in days, and f(p), f(p) are dimensionless func¬ 
tions such that f(p) — f(p) — 1 for p — 2.5 and f(p) = 
10, /(p) = 15 for p — 2.1. All the pre-factors are chosen such 
that they are expected to be close to 1, and therefore the 
leading numerical values are an order of magnitude estimate 
for cb in this regime. It is evident that very low values of cb 


© 2002 RAS, MNRAS 000 , 1 -?? 



6 P. Beniamini et al. 


are required in order for the X-ray and the GeV observations 
to be compatible with each other in the SSC suppressed type 
solution. This is expected, as very weak magnetic fields are 
essential in order to have a sufficiently strong SSC cooling. 


4.2 Slow cooling case 

We assume now that z/ m < vx < v c < ^GeV • In this case the 
X-ray flux is given by (Granot V Sari 2002): 


Fu = < 


0.461 (p — 0.04) e 2,53p (1 + z)^ 


1 +P 


, ,, 2+P 3(1~P) „ 1 -p 

nj E k in ,52 Ws d ~L 2 s for ISM 


5 ±p 


-p- 1 


1 +P 


3.82 (p - 0.18) e 2 ' 54p (1 + z) W e e ; 

t ±p 1 ~ 3 p 0 i-p 

. A * E kin, 52 *days d L28 (j^) 2 mJ y> for wind 

( 6 ) 

Comparing the observed X-ray flux at late times with the 
flux given by Eq. 6 we obtain es(Vi) or e#(A*): 


£b = < 


10 5 g(p) (^) 


2(p — 1) 

(P+4) ~2(p + 3) 2(p+2) 

rp p+4 771 p +4 

r GeV,-6 r X,-5 


~(P+2) (2-3p)(p + 3) 3(p+2) (p— 1) 

P+4 1 2p+8 1 2p+8 

“'0 L GeV,- 3 


d£+ g 4 for ISM 


4 x 10 -t> §(p) ^ 


~(P+ 2 ) 2 — 3p 

A* 1+P 


P + 2 

E?-l pHP 
■ GeV,-6 r X,-5 


-f- A ~h 

b GeV ,~3 L X 


(3p-l)(2+p) 

4(1+P) d7 2 for Wind. 


( 7 ) 

where g(p),g(p) are dimensionless functions such that 
#0) = = 1 for p = 2.5 and y(p) = 4,</(p) = 0.9 for 

p = 2.1. Once again, cb should be very low to allow for this 
type of solution. 


5 NUMERICAL MODELING 

The analytic estimates show that for reasonable external 
density parameters (n ~ 1cm -3 or A* « 1), es should be 
very low in both the SSC suppressed and the slow cooling 
scenarios. Combining the estimates in §4.1, §4.2 we obtain: 
10- 6 < es < 10 -3 . This result is in agreement with recent 
findings (Kumar V Barniol Duran 2009, 2010; Barniol Duran 
V Kumar 2011; Barniol Duran 2014; Santana et al. 2014; ?; 
Wang et al. 2015). 

We turn now to a detailed numerical modeling of the 
afterglow, which might also help us to distinguish between 
the two solutions. We numerically model the afterglow emis¬ 
sion and compare the expected fluxes with the GeV, X-ray 
and optical observations listed in Table 1. In the numerical 
model we relax all the simplifying assumptions made above. 
In particular we relax the assumptions that u c is below the 
GeV band at tcev or that the GeV and X-ray fluxes are 
dominated by synchrotron (and not SSC). 

We calculate the synchrotron and SSC spectra following 
Granot V Sari (2002) and Nakar et al. (2009). In particular, 
the method developed by Nakar et al. (2009) allows us to 
compute the Compton parameter Y(y) as a function of the 
electron energy and to include possible corrections due to 
KN, that may affect the shape of the synchrotron spectrum. 


We ignore the fact that at a given moment in time we ob¬ 
serve emission from different emitting radii (Granot et al. 
1999). We used (Granot et al. 1999; Granot V Sari 2002), 
to verify that the magnitude of this effect cannot change 
the flux as compared with the spectra we are using by more 
than ~ 50%. For each burst, we consider values of e#, n and 
p in the range 10 -8 < cb < 10 -2 , 2.1 < p < 2.8, and 10 -2 
cm -3 < n < 10 2 cm -3 . We keep the fraction of shocked 
energy in the electrons to e e ~ 0.1. We calculate the fluxes 
at tGeV ,tx, and t op t (the latter is the time of optical ob¬ 
servations, in case such observations exist) at the observed 
frequency (100MeV, 1 keV, ~eV respectively). We consider 
a model as acceptable if the difference between the estimated 
fluxes and the data points is within the errors (where the 
errors take into account both the uncertainties of the obser¬ 
vations and of the model). We remark that for some bursts 
it is possible that more detailed observations (i.e. in more 
frequency bands and at different times or using the spectral 
indices in the available observed bands) will further con¬ 
strain the allowed parameter space 2 . However, since we are 
mainly interested in limits, this only means that at worst 
our results may be “too conservative”. 

Optical observations are available for eight out of ten 
bursts in our sample and we include them in the model¬ 
ing. The optical data is missing for GRB 100414A and GRB 
110625A. The X-ray band may either be above or below 
cooling (corresponding to the “SSC suppressed ” and “slow 
cooling” cases respectively). However, there are generally no 
reasonable solutions with the optical band being above the 
cooling frequency at times of order ~ 1 day, with the possi¬ 
ble exception of GRB 090926A for which such a solution is 
only marginally consistent for low values of p ~ 2.1 (see also 
Cenko et al. 2011). For all considered GRBs we can suc¬ 
cessfully reproduce the observed fluxes within the general 
synchrotron or synchrotron self Compton model described 
above. This supports the interpretation of GeV photons as 
radiation from the external forward shock, at least for those 
bursts considered in this study, which are all the bursts with 
a long lasting, power law decaying GeV emission. 

We present, first, detailed results for one GRB and 
then we discuss more generally the results obtained for 
all bursts in our sample. We choose to discuss the case of 
GRB 080916C in an ISM environment. For this burst both 
the “suppressed SSC” and the “slow cooling” solutions can 
account for the observations. In addition, there are some 
solutions at large densities (n > 30 cm -3 ) and very weak 
magnetic fields (es < 10 -6 ) in which both the GeV and 
X-ray fluxes are dominated by SSC instead of synchrotron 


2 In general, errors on spectral indices can be large. Instead of 
using spectral indices, we consider (both for X-rays and optical) 
two data points at times as far as possible one from each other 
(see §5.1 where we do use spectral indices). This assures that we 
are recovering the correct slope of the light-curve. For the GeV 
we consider two data points after T 90 , in order to avoid possible 
contamination from the prompt. We are therefore considering the 
data point where the flux is smaller, and the errors on the spectral 
index and temporal decay can be very large. However, we know 
that LAT spectra that are rising in vF' u are rarely observed and 
in any case their photon index is never harder than -1.7. In order 
to consider only realistic solutions, we add this requirement on 
the LAT spectral index. 
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Figure 2. The allowed parameter space for GRB 080916C in an ISM environment. For any given point in the eg-n plane, colours depict 
the minimum possible value of the isotropic equivalent kinetic energies for which a successful modeling is found (see the values in the 
colour bar). Lower limits on Eq^ u correspond to upper limits on the prompt efficiency (in parenthesis). The left panel corresponds to 
the ”SSC suppressed” solutions discussed in §4 and the right panel to the “slow cooling” solutions discussed in the same section. For 
reasonable values of 3 X 10 -2 cm -3 < n < 3 cm -3 , we obtain: e 7 < 0.55, > 3 X 10 54 ergs and < 5 X 10 5 , independent of the 

type of solution. 




Figure 3. Observed and predicted spectra for GRB 080916C in a homogeneous medium for a single model at six different times: two 
LAT observation times (red, 0.0046 and 0.001 days) two X-ray observation times (blue, 6.94 and 11.57 days) and two optical observation 
times (green, 1.39 and 3.47 days). Left: an example of modelling where the X-ray band is dominated by synchrotron radiation and 
resides below the cooling frequency, while the GeV band is also dominated by synchrotron radiation but resides above the cooling 
frequency. This situation corresponds to a typical “slow cooling” case. The results are shown for eg = 4.5 X 10 -7 , n = 0.14 cm -3 ,^ = 
1.5 X 10 55 erg,p = 2.43. Right: an example of modelling where X-rays are produced by IC-suppressed fast cooling synchrotron emission, 
while the GeV radiation is produced by non-suppressed fast cooling synchrotron emission (corresponding to a typical “SSC suppressed” 
case). Results are shown for cb = 3.5 x 10 _7 ,n = 0.63 cm -3 ,^ = 1.2 x 10 55 erg,p = 2.36. The solid curves are the results of the 
numerical calculation at different observation times: LAT in red, X-rays in blue and optical in green. The observations are denoted by 
filled circles with error-bars (the optical error-bars are smaller than the size of the data points). The red and blue dashed lines are the 
expected synchrotron spectrum above v c , at the time of LAT and XRT observations, under the assumption the flux is not suppressed 
by IC (Eq. 1). It is evident that, according to this modelling, this assumption is correct for the GeV data, while it is not correct for 
X-ray observations, either because the X-ray observations lie below v c (left panel) or because at the X-ray frequency the flux is strongly 
suppressed by IC (right panel). 


emission. However, a more careful examination reveals that 
all the SSC dominated solutions for this burst correspond 
to a GeV photon index T > —1.7 (N u oc v T ) which is not in 
agreement with observations (for all the bursts in our sam¬ 
ple, the photon index is known and it is never harder than 
-1.7. In fact, this is true even when considering the entire 
sample of GRBs detected by GeV, as can be seen in Fig. 


25 of Ackermann et al. 2013a). Fig. 2 depicts the allowed 
parameter space separately for the “SSC suppressed ” (left 
panel) and the “slow cooling” (right panel) solutions. The 
colored region depicts the allowed region in the es-n plane. 
In some cases, for a given pair of values (eg, n) it is possible 
that several solutions are found, corresponding to different 
vales of Eo,kin and/or p. In this case we show the lowest 
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value of Eo t kin, regardless of the value of p. The reported 
values of E 0j kin should then be considered as lower limits. 
The Compton parameter Yx can vary between 3-100. At the 
lower end of this range, Yx is affected by KN suppression 
and may differ somewhat from the expression derived for 
the Thomson regime (Eq. 2) used for the analytic estimates 
presented in §4. Moreover, Yc e v can vary between 0.1-3, im¬ 
plying that a moderate suppression of the GeV flux might 
also take place, but only when 6b is very small and n is 
large. For Yo e v < 0.3 and for X-rays in Thomson regime, 
the numerical results are indeed compatible with the ana¬ 
lytical formulas presented in Eqns. 5, 7. The shape of the 
allowed parameter space can be understood in the following 
way. First, the border to the lower left of the “SSC sup¬ 
pressed ” region is defined by vx — v c at the time of X-ray 
observations. Below this line, the X-ray flux is produced by 
slow cooling electrons (this roughly corresponds to the up¬ 
per right border of the “slow cooling” solutions, where some 
superposition is allowed due to variations of the parameters 
E 0 ,kin and p). The upper border in both cases is due to the 
fact that for larger densities the GeV component becomes 
dominated by SSC instead of synchrotron radiation. Finally 
the upper right border for the “SSC suppressed” solutions 
arises from the requirement that Yx should be large enough 
as compared with Yc e v in order to account for the flux dis¬ 
crepancy between the GeV and the X-rays. All the allowed 
solutions correspond to isotropic equivalent kinetic energies 
satisfying Eo,kin ^ 3x 10 54 ergs, and imply that the prompt 
efficiency is moderate: e 7 < 0.55. This limit on the efficiency 
rules out the value inferred just from the X-ray observa¬ 
tions (e 7 = 0.85) that is derived assuming that X-rays are 
above u c . The allowed parameter space results in strong up¬ 
per limits on the magnetization: < 5 x 10 -5 , regardless 

of the type of solution. Strong limits on the magnetization 
(es < 3 x 10 -4 ) can be obtained also for solutions in a wind 
medium, consistent with previous estimates of the magneti¬ 
zation for this burst (Kumar V Barniol Duran 2009; Gao et 
al. 2009; Zou et al. 2009; Feng V Dai 2010). Typical SEDs 
(spectral energy distributions) for the “SSC suppressed” and 
“slow cooling” cases for this burst (for a particular set of pa¬ 
rameters) are shown in Fig. 3. The solid curves are the SEDs 
resulting from numerical calculation at six different times: 
at the time of the two LAT observations (red), at the time 
of the two X-ray observations (blue), and at the time of the 
two optical observations (green). The six observations are 
denoted by filled circles with error-bars (the optical error- 
bars are smaller than the size of the data points). The red 
and blue dashed lines are the expected synchrotron spec¬ 
trum above is c , at the time of LAT and XRT observations, 
under the assumption that the flux is not suppressed by 
IC (Eq. 1). It is evident that, according to this modelling, 
this assumption is correct for the GeV data, while it is not 
correct for X-ray observations, either because the X-ray ob¬ 
servations lie below v c (left panel) or because at the X-ray 
frequency the flux is strongly suppressed by IC (right panel). 

A similar analysis has been performed for all bursts in 
our sample, for both constant and wind-like density profiles 
of the external medium. For a constant profile (ISM), four 
bursts have both “SSC suppressed” and “slow cooling” type 
solutions, while one has only a “SSC suppressed” solution 
and another has only a “slow cooling” solution. For the re¬ 
maining bursts (090510, 090902B, 100414, 110625) it is im¬ 


possible to find acceptable models to the data, for an ISM. 
For a wind-like medium, no bursts have “SSC suppressed” 
type solutions, whereas seven bursts have “slow cooling” so¬ 
lutions. For the remaining three bursts (090926A, 090902B 
and 090510) a good modeling of the data is possible, but 
it does not correspond to any of the solutions discussed so 
far. In fact, for these bursts the GeV flux is dominated by 
SSC (but still consistent with the requirement on the GeV 
spectral slope: T < —1.7). whereas the X-ray flux is domi¬ 
nated by synchrotron from fast cooling electrons (for GRB 
090926A X-rays may also be dominated by SSC in the solu¬ 
tion for a wind type medium) 3 . 

Fig. 4 depicts the upper limits on es as a function of 
the external density for each burst, for an ISM (left upper 
panel) and for a wind (right upper panel). The bottom pan¬ 
els show the required amplification factor of the magnetic 
field beyond shock compression, denoted AF by Santana et 
al. (2014), for the given cb and n, for a seed field of 10/iG. 
The upper limits are approximately 10 -5 < €b < 10 -3 and 
3 < AF < 1000 for both ISM and wind environments. Lower 
limits on the (isotropic equivalent) kinetic energies as a func¬ 
tion of the external density for both ISM and wind are shown 
in Fig. 5. The limits range between 10 53 — 3 x 10 55 erg and 
are generally much larger than the kinetic energies estimated 
from X-rays only, assuming fast cooling synchrotron 4 . These 
limits on the isotropic energy can be translated into limits on 
the collimation corrected energy if the jet opening angle is 
known (or if a lower limit on the jet opening angle can be es¬ 
timated from the non-detection of a steepening in the after¬ 
glow light curve - i.e., a jet break). The measurements/limits 
on the jet break times are listed in Table 1 (one before last 
column). These have been translated into estimates/limits 
on the jet opening angles (Sari et al. 1999; Chevalier V Li 
2000; Frail et al. 2001). Fig. 6 shows the derived limits on the 
true (collimated) jet kinetic energy as a function of density. 
All GRBs are consistent with an energy < 10 52 ' 5 erg. The 
shapes of the curves in Figs. 4, 5, 6 are quite complex and 
can vary from burst to burst. This is due to the many differ¬ 
ent parameters involved in determining these limits. For low 
enough densities (0.01 cm -3 < n < 1 cm -3 , 0.01 < A* < 1), 
the solutions are usually of the “slow cooling” type. In this 
regime, the upper limit on es is determined by the condi¬ 
tion v c — vx- In some bursts (but for ISM only), for some¬ 
what larger densities (0.1 cm -3 < n < 3 cm -3 ), there are 

3 It is possible that a pure synchrotron solution may be found 
also for GRB 090510 and 090902B; however, it would require 
some changes of the assumptions of the present model, such as a 
more general power-law profile of the external density, a devia¬ 
tion from a single power-law distribution of electrons, etc. Indeed 
for GRB 090510, previous studies proposed synchrotron solutions, 
but these relied on either very low external densities, energy injec¬ 
tion or evolution of the micro-physical parameters (De Pasquale 
et al. 2010; Corsi et al. 2010). For GRB 090902B, Barniol Du¬ 
ran & Kumar (2011) have found synchrotron solutions but with 
a deviation from a single power-law distribution of electrons (and 
even in this case the limits on cb remain very strong). 

4 The only case where the kinetic energy is not much larger than 
that estimated from the X-rays is that of GRB 090328, for an ISM 
and large densities n ^ 1 cm -3 . In this case the solutions are such 
that the GeV is dominated by SSC instead of synchrotron and 
the kinetic energy is similar to that estimated from the X-rays 
assuming fast cooling synchrotron. 
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“SSC suppressed” type solutions, and in this case the limits 
on cb arise from the requirement that Yx should be suffi¬ 
ciently large to account for the ratio of fluxes between GeV 
and X-rays. For larger densities (1 cm -3 < n < 10 cm -3 , 
1 < ^4* < 10) there might be solutions where the GeV flux 
becomes dominated by SSC, and, if such a solutions exist 
they are limited by the fact that as increases the SSC 
peak becomes too weak and the flux in the GeV band is no 
longer sufficiently large as compared with that in X-rays. 
Finally, for the largest densities (10 cm -3 < n < 100 cm -3 , 
10 < A* < 100) the flux may become dominated by SSC 
both in the GeV band and in the X-rays. Since the opti¬ 
cal band in this case is still dominated by synchrotron, the 
allowed solutions are then limited by the fact that as in¬ 
creases the SSC peak becomes weaker and cannot account 
simultaneously for the GeV and optical fluxes. 

In principle, another constraint on these solutions arises 
from the fact that there is no evidence for a passage of the 
cooling break through the X-ray band during the time of 
observation. This constraint is only relevant for bursts in a 
wind environment (for which v c increases with time) and 
“SSC suppressed ” solutions, and for bursts in an ISM envi¬ 
ronment (where v c decreases with time) for the “slow cool¬ 
ing” solutions. However, due to the fact that X-ray observa¬ 
tions are typically available for less than 2 decades in time 
(say 0.1-10 days), and since v c evolves only slowly with time 
(v c oc t -1//2 for an ISM environment and v c oc t 1 ^ 2 in a 
wind), v c cannot change by more than an order of mag¬ 
nitude during the duration of the observation. This is less 
than the spectral range of XRT (0.3-10 keV), and since we 
are using integrated fluxes in this work, it will result in an 
undetectable spectral change given the quality of the data. 

5.1 Results from modelling of simultaneous X-ray 
and GeV data 

For four bursts in our sample, early time X-ray observations 
(simultaneous to observations in the GeV band) are avail¬ 
able. So far we have considered only late time X-ray data, 
because our aim is to test if these data can be safely used 
to infer the blast wave kinetic energy. In this section we ex¬ 
ploit early time simultaneous X-ray/GeV detections (when 
available) in order to test the consistency of our modelling 
and possibly reduce the parameter space identified by the 
modeling of late time X-ray data. 

For all bursts we find that the parameter space overlaps 
with the one inferred from late time observations. In other 
words, it exists a region of the parameter space that allows a 
good description of both early and late time X-ray and GeV 
fluxes. In the following we discuss in detail our findings, for 
each one of these bursts. 

5.1.1 090510 

For this burst we use data at 114 seconds, taken from 
De Pasquale et al. (2010). Similarly to when we use late 
time X-ray data, we do not find a good match between 
the observations and the numerical calculations for a ho¬ 
mogeneous medium. For a wind medium we find values of 
tB = 10“ 7 - 3 x 10" 5 , E 0Mn = 8 x 10 52 - 6 x 10 53 ergs, 
A* =3 — 30. These ranges of values are similar to those ob¬ 
tained with the late time data, but somewhat wider. These 


parameters correspond to a situation in which the X-rays 
are either dominated by fast cooling synchrotron or by SSC 
(for the larger densities) emission, whereas the GeV band 
is always dominated by the SSC component. An example 
of modelling (where the X-ray flux is dominated by syn¬ 
chrotron) is shown in Fig. 7 (upper left panel). The dashed 
lines show the separate contributions from synchrotron and 
SSC, and the solid line is their sum. X-ray and GeV ob¬ 
servations, as well as the errors in the fluxes and spectral 
slopes are denoted by the blue and red error-bars corre¬ 
spondingly. The dotted line is an extrapolation of the ex¬ 
pected synchrotron spectrum above v c (Eq. 1), showing that 
in all cases, the X-rays are incompatible with lying on the 
same segment of the spectra as the GeV photons. These re¬ 
sults are similar to those obtained by Panaitescu (2011), 
who show that the 100 s data is best modelled by a fast 
cooling synchrotron component for the X-rays and an SSC 
component for the GeV. They also find that large isotropic 
energies, large densities and very low cb are required. Other 
studies (De Pasquale et al. 2010; Corsi et al. 2010) have 
shown that the observations can be explained in the con¬ 
text of a homogeneous medium if a very small density is 
invoked (n ~ 10 -6 cm -3 ), well below the range of densities 
considered in this work. These studies also outline that some 
deviation from the standard forward shock scenario (energy 
injection, evolution of the microphysical parameters in time 
or a structured jet model) is required in order to obtain a 
good description of the observations. 

5.1.2 110625A 

For this burst we model the SED at 260 seconds. Both the X- 
ray and GeV data are taken from Tam et al. (2012). Once 
again, as for the late time data, there is no available pa¬ 
rameter space if a homogeneous medium is assumed. For a 
wind medium, however, two sets of parameter ranges are 
found, similar to what was found from the modelling of 
the late time data. The first one (cb = 10 -7 — 6 x 10 -5 , 
Eo,kin = 3 x 10 53 — 3 x 10 54 ergs, A* = 0.6 — 30) corresponds 
to a X-rays either dominated by fast cooling synchrotron 
or by SSC (for the larger densities) and GeV emission al¬ 
ways dominated by the SSC component. The second set of 
parameters is found at lower densities (A* < 0.03), larger 
isotropic energies (E 0j kin = 10 54 - 2 x 10 55 ergs) and similar 
values of cb- In this case, the X-rays are dominated by slow 
cooling synchrotron whereas the GeV range is dominated 
by fast cooling synchrotron. This last type of modelling was 
found also by Tam et al. (2012) who suggest that the cooling 
frequency at t ~ 300 s should be between 10 — 100 MeV. An 
example of such modelling is shown in Fig. 7 (upper right 
panel). 

5.1.3 110731A 

For this burst we use the data at 277 seconds, taken from 
Ackermann et al. (2013b). For a homogeneous medium we 
find t B = 2 x 10“ 6 - 10“ 2 , E 0 ,kin = 10 53 - 2 x 10 55 ergs, 
n < 3 cm -3 . For these ranges of parameters, the X-rays are 
always dominated by synchrotron (either slow cooling for 
lower densities or fast cooling and SSC suppressed for larger 
densities), whereas the GeV emission is dominated either 
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Figure 4. Top panels: upper limits on cb as a function of the external density for ISM (top left) and wind (top right) environments. 
Lower panels: upper limits on the logarithm of the magnetic field amplification factors (AF) beyond shock compression (for ISM in the 
lower left panel and wind in the lower right panel), assuming a seed (unshocked) magnetic field of 10/iG. 
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Figure 5. Lower limits on the isotropic equivalent kinetic energies as a function of the external density for ISM (left panel) and wind 
(right panel) environments. The limits range between 10 53 — 3 X 10 55 erg and are generally much larger than the kinetic energies estimated 
from X-rays only, assuming fast cooling synchrotron. 


by fast cooling synchrotron or by SSC (for larger densities). 
These solutions are similar to those found with late time 
data. For a wind medium we also find allowed parameter 
space but with almost no overlap with the parameter space 
identified with the late time data. This favours a homoge¬ 
neous medium. An example of modelling where the X-ray 
frequency is above v c but the X-ray flux is suppressed by IC 
and where the GeV flux is produced by un-suppressed fast 
cooling synchrotron is shown in Fig. 7 (bottom left panel). 


5.1.4 130427A 

For this burst we use the data at 25,000 seconds. The 
X-ray data is taken from Maselli et al. (2014), and the 
GeV data is taken from Ackermann et al. (2013c). For 
a homogeneous medium we find cb — 10' “ 7 - 2 X 10“ 3 , 
Eo,kin = 3 x 10 53 — 10 56 ergs, n < 2 cm -3 . Here, the X-rays 
are always dominated by synchrotron (either slow cooling 
for lower densities or fast cooling and SSC suppressed for 
larger densities), whereas the GeV emission is dominated 
either by fast cooling synchrotron or by SSC (for larger den¬ 


sities) . The allowed parameter space almost completely over¬ 
laps that obtained with the previously considered data set. 
For a wind medium two sets of parameters are found. The 
first, with cb = 10 -7 — 10 -2 , E 0j kin = 2 x 10 53 — 10 56 ergs, 
A* < 4. These parameters correspond X-rays dominated 
by slow cooling synchrotron and GeV emission dominated 
by fast cooling synchrotron. Another set of parameters is 
found with larger densities (A* > 5), smaller isotropic en¬ 
ergies ( Eo,kin = 10 53 — 5 x 10 53 ergs) and similar values of 
cb- In this case, both the X-rays and the GeV emission are 
dominated by SSC. Both these sets of parameters are found 
also with the previous data set presented above. An example 
of modelling where the X-rays are dominated by slow cool¬ 
ing synchrotron and the GeV by fast cooling synchrotron is 
shown in Fig. 7 (bottom right panel). 


6 CONCLUSIONS 

We have considered all GRBs detected both by LAT and 
XRT, for which the GeV emission lasts longer than the 
prompt and it decays as a power-law in time. We have ana- 
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Figure 6. Lower limits on the true (collimated) kinetic energies as a function of the external density for ISM (left panel) and wind-like 
(right panel) environments. All GRBs are consistent with an energy < IQ 52,5 erg. 





Figure 7. Observed and predicted spectra for the four bursts with simultaneous X-ray and GeV data (for a particular set of parameters). 
In all panels, the dashed lines show the separate contributions from synchrotron and SSC, and the solid line (which is their sum) is 
the total flux predicted by the model. X-ray and GeV observations, as well as the errors in the fluxes and spectral slopes are denoted 
by the blue and red error-bars correspondingly. The dotted line is the synchrotron spectrum above v c as predicted in case negligible 
IC cooling (Eq. 1), showing that in all cases, the X-rays do not satisfy both these requirement, being either below v c or suppressed by 
IC cooling. Top left: SED for 090510 in a wind medium with = 4 x 10 -7 , A* = 8 : E^ in = 3 X 10 54 ergs,p = 2.3 (in this example 
X-rays are dominated by fast cooling synchrotron and the GeV is produced by SSC). Top right: SED for 110625A in a wind medium with 
€b = 1.6 X 10 -6 , A* = 0.02, E^ in = 3.8 X 10 54 ergs,p = 2.1 (X-rays are dominated by slow cooling synchrotron and the GeV by fast cooling 
synchrotron). Bottom left: SED for 110731A in a homogeneous medium with cb = 1.3 X 10 3 , n = 0.02 cm -3 , E^in = 7x 10 53 ergs,p = 2.7 
(X-rays are dominated by SSC suppressed fast cooling synchrotron and the GeV is produced by un-suppressed fast cooling synchrotron). 
Bottom right: SED for 130427A in a homogeneous medium with cb = 2 X 10 _4 ,n = 0.025 cm -3 , E^ n = 1.3 X 10 54 ergs,p = 2.2 (X-rays 
are dominated by slow cooling synchrotron and the GeV by fast cooling synchrotron). 
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lyzed the late time (around 3 x 10 2 s) GeV radiation and the 
late time (around 1 day) X-ray and optical radiation (optical 
observations are available for 8 out of the 10 bursts in our 
sample). Assuming that all these emissions are produced by 
electrons in the shocked circumburst medium we have es¬ 
timated the conditions at this shock and have elaborated 
on the exact radiation process. We modeled the broad band 
observations as radiation from the forward external shock 
under two different assumptions for the external density: 
constant (ISM) and wind-like. We derived synchrotron and 
SSC spectra (accounting for the role of KN suppression) and 
compared the simulated SEDs with observations. While we 
fix e e — 0.1 we have allowed all the other parameters to vary 
over a wide range of values. 

We have found that for all bursts in our sample it is pos¬ 
sible to account for the broad band (optical to GeV or at 
least X-ray to GeV) observations within the forward shock 
model. For each burst we have found a range of possible so¬ 
lutions, which can be classified into two families, depending 
on the position of i/x as compared to v c . In general, low den¬ 
sities and low es (n < 1 cm ~ 3 ,es < 10 -4 ) correspond to a 
high cooling frequency, above the X-ray band. Contrary to 
the electrons emitting X-rays, GeV radiating electrons are 
in the fast cooling regime and vx and VGeV do not lie in the 
same portion of the synchrotron spectrum. In this case the 
GeV flux is a better proxy for the kinetic energy, since the 
X-ray flux depends not only on Eo^kin and e e but also on n 
and cb . Since most of the synchrotron spectrum is produced 
by slow cooling electrons, these solutions require, in general, 
large energies. For larger values of es and/or larger densi¬ 
ties the cooling frequency falls below vx- In this case both 
the X-rays and GeV photons are in the high energy part 
of the spectrum and in order to account for the relatively 
low X-ray flux as compared to the GeV flux, the X-ray flux 
must be suppressed via SSC cooling, while, because of the 
KN suppression, the GeV flux is not similarly suppressed. 
Also in these solutions, the GeV flux is a better proxy for 
the kinetic energy as compared to the inverse Compton sup¬ 
pressed X-ray flux, whose value depends on the Compton 
parameter and hence indirectly on density and es- 

Even if the modeling within the forward shock model is 
successful, only limits to the values of the unknown param¬ 
eters can be inferred, due to the degeneracy among these 
parameters. Wind solutions exist for all the explored range 
of densities (10 -2 < A* < 10 2 ). On the other hand ISM 
requires n < 5 cm -3 . For n and A* in the range 0.1-10, the 
largest possible values for es vary in the range 10 -3 — 10 -5 , 
and the upper limits on the amplification factors of the mag¬ 
netic field beyond shock compression are 3 < AF < 1000. 
The lower limits derived on E 0 ,kin are always larger than 
the estimates derived using the standard assumption that 
the X-ray flux around one day is above v c and that it is not 
affected by Inverse Compton cooling. The inferred values 
of the kinetic energy are sometimes large, but they reflect 
the isotropic equivalent energy, and not the true energy of 
the blast wave. Moreover, we are considering here a sample 
of very energetic bursts, among the most energetic ever de¬ 
tected. Their energetics are not representative of the whole 
sample, since they likely belong to the high-energy tail of the 
population. The lower limits on Eo^in translate into upper 
limits on the efficiency of the prompt mechanism e 7 . These 
upper limits are always less than 50%. 


The values of the magnetic field inferred in our analysis 
are only slightly above the values that would arise from a 
simple shock compression of the magnetic field. It is inter¬ 
esting to note that our results are in line with a recent pro¬ 
posal by Lemoine et al. (2013). In this work Lemoine et al. 
(2013) suggest that shocks amplify the magnetic field signif¬ 
icantly (this is required for particle acceleration within these 
shocks). However, the downstream magnetic field decreases 
rapidly with the distance from the shock. In this case the 
optical and X-ray observations probe values of the magnetic 
field far from the shock. The GeV emitting electrons cool 
rapidly in a region closer to the shock, where the magnetic 
field is still large. However, as the GeV, fast cooling, flux 
is almost independent of the magnetic field, the value of es 
used to model GeV radiation does not affect too much the 
results. A good modeling can be found by using the same, 
small cb for optical, X-ray and GeV observations, with mi¬ 
nor corrections to the GeV flux in the scenario outlined by 
Lemoine et al. (2013). 

The sample considered in this work includes only bursts 
with temporally extended emission detected by LAT. It is 
unclear whether or not these bursts are intrinsically different 
from bursts without GeV emission and represent a separated 
population. However, we note that the bursts detected in the 
GeV range are generally the brightest ones in the sub-MeV 
range (Swenson et al. 2010; Beniamini et al. 2011; Acker- 
mann et al. 2013a). This supports the possibility that this 
kind of emission is typically present in GRBs, but it can 
be detected only for the most fluent bursts, being below the 
LAT detection threshold for fainter GRBs. If this is the case, 
then the bursts studied in this paper are not intrinsically dif¬ 
ferent from the general population of GRBs. If they are part 
of the same population, this triggers the question of whether 
our results (very low e#, X-ray flux either below v c or sup¬ 
pressed by SSC, and prompt efficiencies around < 20%) can 
be extended to the whole population or not. Of course, since 
the requirement of LAT detection selects only the brightest 
bursts, the kinetic energies involved in their afterglow blast 
waves are not representative of the whole population. 

Concerning the tendency to have very low es values, 
we note that there is no reason to imagine a bias towards 
small values introduced by the study of a sample of bright 
bursts. The selection criterion is based on the detection in 
an energy range (~GeV) that lies above ^ c , where the flux 
is very weakly dependent on es’- F v oc e^ -2 ^ 2 . This implies 
that if a selection bias exists, it goes in the opposite direc¬ 
tion, since the GeV afterglow flux would be larger for higher 
values of es- Moreover, two GRBs with the same kinetic 
energy but very different cb would have very similar GeV 
fluxes, and then the same probability to be detected by LAT. 
We also note that small values of es have been derived also 
in recent studies of larger samples of bursts without GeV 
radiation (Barniol Duran 2014; Santana et al. 2014; Zhang 
et al. 2015; Wang et al. 2015; see also the exercise of not 
including GeV data in the afterglow modelling of bursts de¬ 
tected in the GeV in Kumar V Barniol Duran 2009, 2010; 
Barniol Duran V Kumar 2011; Lemoine et al. 2013). It is 
not easy to understand why recent works (with and without 
LAT data) are highlighting the existence of cb values much 
smaller that those derived in previous analysis. We can spec¬ 
ulate that, at least in cases where data are not sufficient to 
break the degeneracy between the several free parameters 
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(especially cb and n) solutions corresponding to very small 
cb values are usually discarded, in favour of values consid¬ 
ered more standard. Moreover, also in those cases where the 
free parameters are all constrained by the data, addition of 
GeV data obviously brings some new information and opens 
a region of the parameter space which is usually not investi¬ 
gated when high-energy observation are missing. In fact, as 
we have proved in this work, GeV data are inconsistent with 
being a simple extrapolation of the X-ray data, and this rules 
out large cb values, in favour of smaller values. Small values 
of es require a proper treatment of the IC cooling and KN 
effects. These effects modify the shape of the synchrotron 
spectrum and allow us to recover a satisfactorily modelling 
of the broad band data and to explain the apparent inconsis¬ 
tency between X-ray and GeV fluxes. Our results, together 
with results from samples without LAT detection, suggest 
that the distribution of es in the forward external shocks 
of GRBs is wider than what is usually believed, possibly 
rangingall the way down to 10 -7 . 

The availability of LAT data affects not only the esti¬ 
mate of €b , but also the estimates of the blast wave kinetic 
energy (and then, the estimates of the prompt efficiency). 
Also in this case, it is interesting to compare estimates de¬ 
rived by including and excluding LAT data. In our sample, 
the efficiencies derived from X-ray alone agree with results 
found for other bursts, i.e. large values close to unity (§3). 
Again, only the inclusion of LAT data forces the modelling 
to prefer higher values of the kinetic energy, lowering the 
efficiency requirement on the prompt mechanism. The fact 
that the GeV detected bursts, which are most energetic in¬ 
dicate low efficiency during the prompt phase, suggests that 
the results are generic. Otherwise why would the brightest 
burst be also the least efficient? 

Regardless of whether or not these results can be ex¬ 
tended to a large sample of bursts, our study outlines the 
need to perform broad band afterglow modeling in order to 
properly infer important properties like the energetics and 
the efficiency, and put constraints on the nature of the emis¬ 
sion mechanisms in GRBs. 
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7c = l + y = ^ h „<7tB 2 (1 + Y) = 4w(l + y)e B A J R'<7 T ’ 

(A.4) 

where u — rn' 0 m p c 2 is the energy density of the shocked 
fluid, t' dyn is the time since the shock began and A R' is 
the shell’s width. The Compton Y parameter in Eq. A.4 is 
actually Y c , which corresponds to electrons radiating at v c \ 
however, for all electrons in the Thomson regime, Y is the 
same. 

For 7 C > 7 m , Eq.A.2 yields n' 0 « Cjm/{p - 1) and we 
can write the average Lorentz factor squared as: 



dn' e 2 , 

— 7ed7e = 
d^fe 


(p - l)7m 7c P _ 
(3-p)(p-2) 

(A.5) 


The Compton parameter is given by: 



4Q - 1 ) 777 ^ x y^ p _ 
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Plugging in r e = nocrrAii', and using Eqs. A.3, A.4 we 
obtain: 
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Notice that in this equation, Y appears also indirectly 
through v c (see Eq. 4). 


APPENDIX A: SSC COMPTON PARAMETER 
IN THE THOMSON REGIME 


We give a short derivation of Eq. 2 used for estimating the 
SSC Compton parameter in the Thomson regime. We as¬ 
sume 2 < p < 3 and Vm < which is the most relevant 
situation for GRB afterglows at times between 10 3 — 10 5 sec, 
when late GeV and late X-ray / optical observations are 
taken. Assuming an initial power law spectrum of electrons 
with a slope p , and taking into account evolution of the par¬ 
ticle spectrum due to cooling we write the density per unit 
energy as: 

dn' e _ ( C(^e/jm)~ P , 7m < 7e < 7 c 

dye - \ C'(7c/7m) _p (7e/7c)“ p “ 1 , 7c<7e 

(A.l) 

where here and throughout this section, primes denote quan¬ 
tities in the comoving frame. C is given by the total density 
of electrons, n ' 0 : 
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Next, we relate 7 m ,7c to e e ,€B'- 
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